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SUMMARY 


An analytic investigation is made of the aerodynamic coefficients of 
various todies located in a free-molecule-flow field. These todies are 
the following: flat plate, cylinder, sphere, and cone. Calculations are 

performed using values of molecular speed ratio (ratio of stream speed to 
most protatle molecular speed) ranging from 0 to 20. 

* The aerodynamic coefficients of a cone are calculated for angles of 
attack ranging from 0° to 60°. The semi vertex angles of the cones investi- 
gated vary from 2.5° to 30 °. 

* 

The calculations are performed assuming two typeB of molecular 
reflection, specular and diffuses in addition, for the cone, a third type 
of molecular reflection, wherein impinging molecules are not re— emitted 
from the body hut are swept along its surface, is postulated in order to 
compare the drag coefficients calculated by free-molecule— f low theory with 
values obtained from continuum theory. 


INTRODUCTION 


Consideration of the problem of the flight of high-speed long-range 
aircraft has Indicated that skin temperature and drag forces may be 
considerably reduced by operation at high altitudes. Consequently, it 
becomes of interest to determine the aerodynamic forces arising from high- 
speed flight at high altitudes. Flight at altitudes where the molecular- 
mean— free path is larger than a characteristic body dimension is in the 
free-molecule-f low regime . Molecular-mean— free path as a function of 
altitude is given in reference 1. (At an altitude of 75 miles, the mean- 
free path is about 1 foot.) Calculations relating skin temperature, 
altitude, and velocity have previously been completed in reference lj 
drag forces, however, were not computed. 

Analytic investigations of the drag forces acting on simple body 
shapes in free-molecule flow have been reported in references 2, 3* and 
4, while experimental, investigations have been described in references 
4 and 5« 
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In all these analyses, it was assumed that the gas molecules have a 
Maxwellian distribution of thermal velocity superimposed upon the mass 
velocity. This assumption is likewise made in the present paper. 

Beineman (reference 2) and Ashley (reference 3) bave calculated the aero- 
dynamic coefficients for various body shapes including the flat plate, 
sphere, cylinder, and cone. However, the expressions for the drag coef- 
ficients of the sphere and cylinder were not integrated by Heineman, 
while Ashley’s solution for the cylinder was not expressed in closed 
form. Also, in both papers, the drag coefficients for the conical bodies 
were calculated only for the special case of zero angle of attack. In 
this paper, the drag coefficients for the sphere and cylinder are ex- 
pressed in closed form, and the calculations for the drag coefficients 
for the cones are extended to include angles of attack other than zero. 
Two types of molecular reflection, diffuse and specular, are investi- 
gated. Molecular reflection which occurs in a random direction without 
relation to the previous velocity direction is termed diffuse reflec- 
tion. Reflection such that the molecules leave the surface at an angle 
equal to the angle of incidence is called specular reflection. 

As previously stated, free-molecule theory accounts for both the 
mass motion of the gas and the thermal motion of the molecules. Early 
investigations of rarefied gas flow were based upon the assumption that 
the velocity of the molecules was equal to the translatory velocity of 
the gas, thus neglecting the thermal motion of the molecules. As 
pointed out by Zahm (reference 6), expressions for the forces acting on 
various bodies located in a field of such particles were derived by 
Newton, who assumed that the molecules were either perfectly elastic or 
perfectly inelastic. A perfectly elastic molecule rebounds after a 
collision with a surface in such a manner that the n rrnal component of 
velocity is reversed. A perfectly Inelastic molecule loses the normal 
component of velocity upon impact. In both, cases, however, the tangen- 
tial component remains unchanged. In this paper, flow of elastic par- 
ticles having only a translatory mass motion will be termed "Newtonian 
flow," and flow of inelastic particles will be termed "Inelastic 
Newtonian flow." 

At high flow speeds where molecular thermal velocities become rela- 
tively unimportant, the expressions derived from free-molecule— flow 
theory for the forces on a body for the case of specular reflection re- 
duce to the expressions for a body in a Newtonian flow field. 

7 .'.elastic Newtonian flow approximates the condition of continuum 
hypersonic flow in which the shock wave lies very close to the body. 
Consequently, impact forces of inelastic Newtonian flow have been applied 
to an analysis of hypersonic flow (reference 7). In addition to the 
impact forces of Inelastic Newtonian flow, other forces of a centripetal 
nature exist as a result of flow over curved surfaces. An Investigation 
of the effects of these centripetal forces is also presented in refer- 
ence 7« 
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The purpose of this paper is to obtain solutions for the aerody- 
namic coefficients of various bodies in a fre e-mcle cule— flow field and 
to calculate the lift and drag coefficients for a conical body at an 
angle of attack, assuming two types of molecular reflection, specular 
and diffuse. A comparison is then made of the drag coefficient for a 
cone at zero angle of attack calculated by free-molecule methods, assum- 
ing a hypothetical type of molecular reflection similar to that of 
Newton's inelastic particles, with the drag coefficient tabulated in 
reference 7 for inelastic Newtonian flow and the drag coefficient calcu- 
lated in reference 8 for continuum flow. 


NOTATION 


A 

* 

c x 

c z 

- 

Ci 

C Di 

°Li 

C r 

C Dr 



C 


D 


C L 



surface area, square feet 

components of total velocity, feet per second 

aerodynamic coefficient due to impinging molecules, dimensionless 

drag coefficient due to impinging molecules, dimensionless 

lift coefficient due to impinging molecules, dimensionless 

aerodynamic coefficient due to diffuse re-emission of the 
molecules, dimensionless 

drag coefficient due to diffuse re-emission of the molecules, 
dinens lonle s s 

lift coefficient due to diffuse re-emission of the molecules, 
dimens lonle s s 

total drag coefficient assuming diffuse re— emission 
(C Dl + Cj^), dimensionless 
total lift coefficient assuming diffuse re-emission 
(°Li _r dimensionless 

total drag coefficient assuming specular reflection of the 
molecules, dimensionless 
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s r 




total lift coefficient assuming specular reflection of the 
molecules, dimensionless 

momentum force due to impinging molecules, pounds 
normal diffuse re-emission momentum force, pounds 
total normal specular momentum force, pounds 
modified Bessel function of first kind and zero order 
modified Bessel function of first kind and first order 


direction cosines in an arbitrary direction, dimensionless 
(Subscripts l and d refer to lift and drag directions.) 


length, feet 

Mach number, dimensionless 

mass of one molecule, slugs 

number of molecules per unit volume of gas 

radius of body, feet 

molecular speed ratio (ratio of stream mass velocity to most 
probable molecular speed) 


( — = / 2- M ), dimensionless 

W / 


diffuse re-emission speed ratio 
dimensionless 


f u\ 

Iri’ 

\ v r/ 


mass velocity, feet per second 


components of mass velocity, feet per second 
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Tr 


x 

7 

z 


1 


a 


0 


7 

6 

0 

P 

9 

X 


X* 


ft 


ft* 


most probable molecular speed, feet per second 
most probable re-emission speed, feet per second 

local Cartesian coordinates 


angle of attack of "body with respect to free stream 
reciprocal of most probable molecular speed 
seconds per foot 


f 1 \ 

W’ 


ratio of specific heats, dimensionless 
semivertex angle of cone 

angle of attack of body element with respect to free stream, 
density of gas stream, slugs per cubic foot 
azimuthal angle measured at base of cone 
dimensionless quantity defined by 
X = exp (— a 2 Z X( ^ 2 ) + 8 ^xd Cl + erf (s Z^a)] 

dimensionless quantity defined by 

X* = exp C-s 2 ^’ 2 ) - \/rt~ s Cl - erf (s Zxd* )] 

dimensionless quantity defined by 

ft = U z pU 2 + lyPUy + l z f3U z ) |exp(-p £ U s 2 ) + 

a/ 5T PU X [1 + erf(pu x )]| + -5^1 Zx [1 + erf (pU x )] 
dimenBionleBS q uan tity defined by 
ft* = C Z z * 3Z7 z * + Zy’pITy* + Z z ’{3U z ') |e3p(-p 2 U x * 2 ) - 

gUx* Cl - erf (BU r * )]j- — ^Z x * Cl - erf (pUx* )] 
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©rf (a) 


error function 





exp(— x®) dx 


Superscript 


refers to shielded surfaoe only 


ANALYSIS 


The hasic assumptions involved in the formulation of the expres- 
sions for the force on a tody in a free-molecule— flow field are the fol- 
lowing: (1) the gas molecules have a Maxwellian velocity distribution 
superimposed upon a uniform mass velocity, and (2) collisions between 
impinging and re-emitted molecules are negligible. 

As a result of the second assumption, the total force can be broken 
down into two components: one arising from the bombardment by impinging 
molecules and the other arising from the re-emission of the molecules 
from the surface. 

The type of re-emission of the molecules from the body surface can 
be divided into several categories, the more common of these being dif- 
fuse and specular reflection. The phenomenon of molecular reflection Is 
•treated in greater detail in references 1 and 9, and only a brief dis- 
cussion Is given here. Diffuse reflection, which is most common physi- 
cally, occurs in such a manner that all previous directional history Is 
erased, the actual direction of re-emission being controlled by the 
Khudsen cosine law. The nature of this type of re-emission is such that 
the molecules, upon leaving the surface , have a Maxwellian distribution 
of speed which depends upon the temperature of the re-emitted stream. 

For the case of specular reflection, the molecules leave the surface In 
a direction determined by the angle of incidence. Thus, the normal corv- 
ponent of velocity is reversed, while the tangential component remains 
unchanged. 

The method of calculation of and C r , the incident and re- 
emission aerody nami c coefficients, is outlined in reference 4, and Is 
summarized in the following paragraphs. 

In the integrals presented below, which are taken from reference 4, 
t he direction cosines between the local coordinate axes of the body sur- 
face area and the momentum force are given 2y, and lz» Figure 

1 lllustiates the coordinate axes for an element of area used In this 
report. The positive x axis is normal to the element and is directed 
into the exposed surface. The mass velocity vector mak es an angle 0 
with tlie positive direction of the y axis. Figures 2 and 3 illustrate 
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the elements of area chosen for each of the "bodies. The components of 
mass Telocity TJ X ., Tly, and Uz are in the direction of the local axes. 

The first integral term applies directly to the area exposed to the mass 
Telocity; the second integral term applies to the rear or shielded sur- 
face of the "body. These areas are taken into account separately "because 
the front of the "body acts as a shield preTenting collisions of the rear 
side with al 3 molecules except those haTing absolute Telocity components 
in the direction of flight with magnitude equal to or greater than the 
flight speed (reference 1). The shielded area on the surface of a cone 
is shown in figure Other methods haTe "been deTised to account for the 
shielded area and are presented in references 2 and 3* 


Force Due to Impinging Molecules 


The component of momentum in a direction defined "by Z Zj lj, anil 
l z imparted to a differential plane area in a fre e— molecule— i flow field 
"by the impact of impinging molecules is 


cLGl = 

1 it 3 / 2 


/-ico oco 000 r r 

J J J Cx(^x°x + V°y + ^z c z) exp -j— fl 2 (c x — U x ) 2 + 


(cy — Uy) 2 + (c z — U Z ) 2 j" dC X dCy dC Z — 


poo pco no C r 

J J J + ^y°y + T'T.cz) exp — P 2 (ex ~ Ux)‘ 


(cy — Uy) 2 + (c z — U Z ) 2 J j- dc z dCy dc s 


dA 


( 1 ) 


or 
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dGi . fs! i 


( 2 z pU z + 2 yPU y + 2 z 3 IT z ) 0U X [1 + erf (0U X )] + 


2 v^t B ; 

©xp (h 3 2 U z 2 )} +^| 2 X [1 + erf 


>{v^ 


( 2 z f3U x + Z y £U y + Z z 3U Z ) /,/? pu x [l - erf (j3U z )3 - exp (~£ S U X Z ) j- + 


Z x [ 1 - erf (PU X )] 


dA 


( 2 ) 


where P = mR and £ = l/V m = s/CJ. 

The choice of the direction cosines is related to the component 
of force "being computed. If the drag is "being calculated, it is nec- 
essary to find the direction angles between each of the local axes 
and the direction of drag. Likewise, direction angles are found for 
any other farce component "being investigated. 

In applying equation (2) to arbitrary bodies, the momentum force 
in a desired direction due to the impinging molecules is first found 
for an element of area of the body under investigation. "When this 
relationship is used, an integration over the surface area of the body 
gives the desired momentum force on the entire body due to impinging 
molecules. The incident aerodynamic coefficient referred to a char- 
acteristic body area is then found. 


Force Due to Diffuse Molecular Reflection 


For the case of diffuse molecular reflection, the normal 
re-emission force on the front side of a differential area inclined 
to the direction of flow Is given by 


dG r = dA 

r 2 2SST 


and the normal re— emission force on the rear side is given by 


( 3 ) 


dG r » 


pU g X' 

2 2ss r 


w 


dA 
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The difference obtained by subtraction of equation (4) from equation ( 3 ) 
is the total re— emission force. (See appendix A, reference 4.) 

The value of s r is a complicated function of the surface condi- 
tion of the body; however* in this paper, s r is assumed equal to s 
for all calculations. At low flow speeds, this assumption is very 
nearly correct. At high flow speeds, however, the value of s r is less 
than that of s. 

The resolution of the normal force of equations ( 3 ) and (^) into 
the lift and drag components followed by integration over the body sur- 
face gives the re-emission lift and drag acting on the body. After 
calculation of the lift and drag, the aerodynamic coefficients for dif- 
fuse molecular reflection can be found. 


Force Due to Specular Molecular Reflection 


j* For the case of specular reflection, the molecules leave the sur- 

face in a direction determined by the angle of incidence. Since the 
reflection angle is equal to the angle of incidence, the normal component 
r of velocity is reversed, while the tangential component remains un- 

changed. Therefore, the momentum change is in a direction normal to the 
bombarded surface, and the total momentum force consists of two equal 
components; one due to Impingement of the molecules, the other due to 
the reflection of the molecules from the surface. 

The incident momentum force on a differential area in a direction 
normal to the surface is found from equation ( 2 ) by the substitution of 
the following direction cosines; 

Z x = cos 0=1. 

2 y = cos — = 0 ( 5 ) 

, it « 

Ir, = cos — = 0 

z 2 

The total normal momentum force for the case of specular reflection is 
twice as great as the incident normal force, and it is given by 

dGg = ^ |v^T[l + erf (j3U z )] (Vu ^ 2 + |) + )} - 

[1 - erf(3U x )3 (Vu x 2 + “ ^sp (-g 2 U s 2 )| cLA (6) 
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Integration of equation (6) over the surface of the "body after the nor- 
mal force has been resolved into the lift and drag components yields the 
total lift and drag on the body due to specular reflection. 


Force Due to Molecular Reflection Analogous to Hypersonic Flow 


A third type of molecular reflection is postulated here in order to 
approximate the characteristics of hypersonic flow over a body. This 
type of reflection is much like that of Newton's inelastic partioles, 
the only difference being that thermal velocities are considered. It is 
assumed that the molecules flow along the body after collision with the 
surface; thus, the normal component of velocity is destroyed while the 
tangential component remains unchanged. The values of the aerodynamic 
coefficients thus obtained are one— half the values calculated for specu- 
lar reflection. This type of molecular reflection is extremely unl ikely 
and probably physically impossible; however, the Justification for this 
assumption lies in the fact that the conditions of hypersonic continuum 
flow over a body are approximated. In continuum flow, the stream fol- 
lows the surface of the body. As the flow speed increases, the zone of 
influence of the body on the gas stream approaches the surface of the 
body until the shock wave lies nearly on the body surface. On the basis 
of the above assumption, free-molecule theory is compared with continuum 
theory. 

The expressions for the aerodynamic coefficients of a flat plate, 
cylinder, sphere, and cone together with the details of the development 
are presented in the appendix. 


DISCUSSION AND RESULTS 
Flat Plate 


The aerodynamic coefficients for a flat plate inclined at various 
angles of attack to the direction of mass flow were calculated for 
values of molecular speed ratio varying from 0 to 20 and are presented 
in figures 5, 6, and 7, 

The limiting values, as the molecular speed ratio approaches infin- 
ity, of -yie aerodynamic coefficients for a flat plate for the case of 
specular reflection reduce to the classical expressions for Newtonian 
flow, mentioned by flabm (reference 6) 

Cpg = 4 sin 3 a 

C^g = 4 sin 2 a cos a 
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A comparison of the drag coefficients for diffuse an d specular 
reflection indicates that the value for diffuse reflection has a larger 
value at low angles of attack; at the high angles of attack, the spec- 
ular drag coefficient is greater. This is explained by the fact that, 
at small angles of attack, the drag for the case of specular reflection 
is due to the snw.1 1 normal component of velocity; whereas the drag for 
the case of diffuse reflection has a value which results from the loss 
of "both a small normal velocity component and a large tangential velo- 
city. At large angles of attack, the re-emission velocity arising from, 
diffuse reflection is less than the re-emission velocity due to specular 
reflection, while the momentum effects of the incident molecules are 
nearly the same for "both cases. Thus, the total drag coefficient due to 
specular reflection is greater. This last point is made clear "by consid- 
ering a flat plate at a 90° angle of attack. 

The lift— drag ratios for the plate are shown in figure 6. The ex- 
pression for the lift— drag ratio far specular reflection is cot a. 
Obviously, these values have a great range. The aerodynamic coefficients 
for double-wedge wings can easily he found by modifying the coefficients 
for the flat plate. The low values for diffuse reflection are notable. 


Cylinder 


# 

Values of the total drag coefficient for the case of diffuse re- 
flection were obtained from, reference 4 and are presented in figure 8. 
An analytical expression for s r employed in the calculations of ref- 
erence 4 was found by Beans of an energy and mass balance over the sur- 
face of the cylinder taking into account conditions which exist at the 
surface of the cylinder. Details of the calculations can be found in 
appendix A, reference 4 . 

A recent experimental investigation of drag (reference 4 ) has sup- 
plied the drag-coefficient data plotted in figure 8. Theoretical 
curves of the drag coefficient for both diffuse and specular reflection 
are plotted for comparison with the experimental data. A complete dis- 
cussion of the comparison of the experimental data and the curve repre- 
senting diffuse reflection can be found in reference 4 . In addition, a 
curve of the drag coefficient for the case of diffuse reflection, which 
includes the assumption that the molecular speed after impact is the 
same as that before impact, is shown for comparison with the drag co- 
efficient obtained from reference 4 . The curves agree well at the low 
speeds, but diverge as the molecular speed ratio increases. 
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Sphere 


The results of the calculations for the drag coefficient for a 
sphere are presented in figure 9» The drag coefficient for the case of 
diffuse reflection is compared with the drag coefficient for the case of 
specular reflection. The value of the drag coefficient for diffuse 
reflection lies above the value for specular reflection but both curves 
asymptotically approach the limiting value of 2 at large values of s. 

The investigation into the case of specular reflection for a 
sphere yields an interesting result. Equations (A15) and (A19) indi- 
cate that the total drag due to the specular type of reflection has the 
same value as the drag caused by molecules striking the surface and 
coming to rest; that is, the integrated effect of the re-emission drag 
arising from specular reflection of the molecules from the surface is 
zero. It should be noted that the sphere is unique In that the re- 
emission momentum force of the molecules reflected In a forward direc- 
tion is equal and opposite to the momentum force of the molecules re- 
flected in the direction of the mass flow. 


Cones 


Calculations for the drag and lift coefficients of the various 
cones, assuming three types of molecular reflection, are shown in 
figures 10, 11, 12, 13 , and l4. 

Results considering diffuse reflection .— Figure 10 gives the drag 
coefficient, Cp = Cp^ + Cp^, for the case of diffuse molecular re flee— 

tlon. It is interesting to note that the drag coefficients for the 
sharp cones are higher than those for the blunt cones. The impinging 
molecules lose both the tangential and normal components of velocity. 

The tangential or friction force on the greater exposed surface area of 
the sharp cones Is greater than the tangential force on the blunt cones. 
Hence, the drag coefficients are large for the Bharp cones since the 
reference area Is the same base area far el 3 cones. 

The lift coefficient, Cp, = Cp^ + Cp^, Is plotted in figure 11. 

The lift coefficient for blunt cones at large angles of attack de- 
creases as molecular speed ratio becomes small. This is due to the col- 
lision of the top or shielded surface of the cone with molecules posses- 
sing an absolute velocity in the direction of flight with magnitude 
equal to or greater than the flight speed. The momentum force due to 
these molecules is then in a direction opposite to the direction of 
lift. 
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At low values of molecular speed ratio, in tie range -where the 
ther mal velocity is equal to or greater than the mass velocity, the 
contribution of the thermal velocity to the absolute velocity becomes 
a determining factor in the calculation of the aerodynamic coeffi- 
cients j whereas, at high values of molecular speed ratio, the contri— ' 
bution of the thermal velocity is negligible. 

Besults considering specular reflection .— Figures 12 and 13 give 
the drag and lift coefficients for cones for the case of specular re- 
flection. At smal 1 angles of attack, the drag coefficients of the blunt 
cones have higher values than those of the sharp cones. The reverse is 
true at large angles of attack. A very sharp cone at zero angle of 
attack has very little drag because the change in momentum of the mole- 
cules has only a very small component in the direction of drag. ThiB 
component, acting on a blunt cone, is much greater. At a high angle of 
attack, the surface area exposed to the flow by the sharp cones exceeds 
-the exposed surface area of the blunt cones, for equal base areas, and 
hence the drag force is greater for the sharp cones. Since the drag 
coefficients are referred to the same base area rather than a projection 
of the exposed area, they must vary in direct proportion to the drag 
force . 


Besults considering an approximation of continuum theory .— The 
values of the lift and drag coefficients obtained by means of the hypo- 
thetical type of re-emission postulated in the analysis are half the 
values obtained for the case of purely specular reflection in which the 
normal component of velocity is reversed as a molecule leaves the sur- 
face . 


Values of drag coefficients for cones in continuum flow have been 
calculated by Kopal (reference 8). Likewise, reference 7 describes an 
investigation of continuum flow at hypersonic speed. In figure 1*4-, the 
drag coefficient of a 30 ° semi vertex angle ccne, at zero angle of 
attack, calculated for the above hypothetical case of molecular re- 
flection, is compared with the drag coefficient calculated in references 
7 and 8. ■■■— 

Experimental values for the drag coefficient of a 30° semivertex 
cone-cylinder fired in the Ames supersonic free— flight tunnel (reference 
10) are also shown in figure 1*4- for comparison with the ana l ytic 
values. The experimental values include head drag on the cone, friction 
drag on the cylindrical afterbody, and base drag. The free-molecule 
values take into account head drag on a 30 ° semivertex angle cone only. 
At high values of flow speed, the drag coefficient values obtained from 
free-molecule theory approach the values for inelas J Ic Newtonian flow. 
Although the analysis of reference 7 for flow past a cone takes into 
account centrifugal forces arising from flow over curved surfaces, the 
expression for the drag coefficient reduces to that for inelastic New- 
tonian flow at zero angle of attack where the centrifugal forces are 
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zero. The drag coefficient due to inelastic Newtonian flow is repre- 
sented in figure 14 by a straight line., since it is independent of flow 
speed. 

In figure llj-, the agreement between fre e-molecule— flow theory and 
the experimental data for a continuum flow at low values of molecular 
speed ratio is entirely fortuitous. The increase in the value of the 
drag coefficient calculated from free-molecule— flow theory is due to the 
fact that the drag force is proportional to the first power of the velo- 
city at low speeds; whereas the increase in the value of the experimental 
drag coefficient is the characteristic rise which occurs in the tran- 
sonic range of stream speeds in continuum flow. 


CONCLUDING REMAKES 


The lift-drag ratios of a flat plate in free-molecule flow are 
very small for the case of diffuse molecular reflection. The lift- 
drag ratios decrease as the molecular speed ratio is increased. At 
high speeds, the lift-drag ratios for the flat plate are more favorable 
at small angles of attack than at large angles of attack. 

For the case of diffuse molecular reflection, the lift and drag 
coefficients of all the bodies investigated approach constant limiting 
values at high speeds. These limiting values depend upon the body 
geometry and the angle of attack. At low speeds, the drag coefficients 
for all the bodies approach an infinite value; this is due to the fact 
that the drag force is proportional to the first power of the velocity 
at low speeds. 

For the case of specular reflection, the aerodynamic coefficients 
of all the bodies investigated likewise approach constant limiting 
values at high speeds. These limiting values correspond to those cal- 
culated by means of inelastic Newtonian flow theory. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., May 9, 1951. 
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APPENDIX 

DERIVATION OF THE AERODYNAMIC COEFFICIENTS FOR 
BODIES IN A FREE-MOLECULE-FLOW FIELD 


Tlie details of the development for the lift and drag coefficients 
for the various bodies are given below: 


Flat Plate 


The combined form of equation (2) is. 
dGi = {2 (ixWx+ljWj+lz&Uz) v^T pU x erf (pUxJ+expC-B^Ux 2 ) 


Jit S 

Jk Z x erf(PU x )j“ dA 




(Al) 


Since the faces of a flat plate are parallel, the same direction 
cosines apply to both the front and rear surfaces and the momentum of 
the molecules striking both faces can be collectively taken Into account, 
and equation 2 written in the combined form is applicable. This Is not 
true for an element of area chosen for the cone* however, and the com- 
bined form of equation (2) cannot be used for that case. For the case 
of the flat plate, equation (Al) gives the incident momentum force on 
both sides of the plate, and the aerodynamic coefficients can be found 
without further integration. 


The direction cosines and velocity components are the following: 
Z x <i = sin a Z x ^ = cos a U x = U sin a, 

Ijl = sin a TJy = -D cos a 

Izl = 0 Uz = 0 

PD = s 


Z — —cos ct 
^zd = 0 


(A2) 


Substitution of the relations of eq.ua ti on (A2) into equation (Al) 
gives the following coefficients, referred to the area of one side of 
the plate. 
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/■* 2 2 

C;d = — = — exp (-s sin a)+2 sin a 

1 </* B 

cos a 

Cl* = S- srf - ( s 

■*■ s 

The normal force due to diffuse molecular re— emission from both 
front and rear surfaces is given by the difference between equations (3) 
and ( 4 ). 

dGr - dGr‘ = 555- (X -X') dA (A 5 ) 


( 1 + 5pO 


erf (s sin a) (A 3 ) 


sin a) 


(A 4 ) 


Substitution of U sin a for Ux and s for PIT in X and X* 
reduces equation (A 5 ) to 

pU 2 J* sin a . 

dG*. - dGr* = dA (A6) 


The re— emission drag force on an element of area is 

pU 2 4/ir sin 2 a 
~2 


dA 


and 


Cl* - 


sin 2 a 


s r 


(A 7 ) 


The total drag coefficient is then 

o , i . . . . JT sin 2 a 

Cp = exp (~s 2 sin 2 a) + 2 sin a (l + — -g) erf (s sin a)+ — — 

Jn e 


The lift due to molecular re-emission on an element of area is 
pU 2 J~x sin a cos a 

2 s r 


dA 


and 


= 


sin a cos a 


(A8) 


The total lift coefficient is then 


cl 


cos a 
"IS- 


erf (s sin a) + 


nfri sin a cos a 
s r 
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For the case of specular reflection, the total momentum force in a 
direction normal to the surface resulting from inoident and reflected 
molecules on loth sides of a flat plate is 


dG s 


ptf l 2 2 

2 JIT s 


2 £U x exp( -£ 2 TJ Z 2 ) + 




from equation ( 6 ). 


erfOUx) 


dA 

(A9) 


The drag and lift coefficients are then found after resolving the 
normal force as given in equation (A 9 ) into drag and lift components. 

if 

sin a, exp(~s 2 sin 2 a) + 


Cpg = sin a 


4 ^sin 2 *! 


+ ) erf(s Sin a) 

2 s 2 / 


(AID) 


C^g = cos a 


sin a exp(— e 2 sin 2 a,) + 


L b 

4 (sin 2 a + 3*) SrfC s sin a) 


(All) 


Cylinder 


The case of diffuse molecular reflection on a transverse cylinder 
has "been investigated in reference k, and the total drag coefficient 
based on the projected area of the cylinder is 


•J* f B Z \ J /B 2 \ 

= — exp [--zj |I 0 (-z-J 


J D 


(l+ 2 s 2 ) 


(?)*“(?)]} 


US / 2 

^Sv. 


The specular case only remains tj he considered here. The element 
of area is chosen as shown in figure 2, thus 

dA = L E d 0 


l xd = sln 9 

lyi = — 00 s 6 


1 J X = U sin G 
U y = -Or cos 6 


*zd = 0 


U z = 0 
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Each element of area is considered to hare a rear or shielded area 
and the normal force on the cylinder is found from equation (A9). 


dGs = 


pu 2 km 


S 2 ■ L- 


2s sin 0 exp(— s 2 sin 2 0) + 


2 Jn ^ 


s 2 sin 2 S + — ) erf(s sin 0) d0 


The drag force is 

w pu 2 4lr rf 

2 JTL s 2 J 0 


sin 0 2s sin 0 exp(-s 2 sin 2 0) + 


2,v£ ^s 2 sin 2 0 + erf(s sin^0) d0 (A12; 

The integrals are 

sin20 exp(-s 2 sin 2 0) d0 = | exp ^ jlo ( jp) “Ii (■'§“) J 

jr/2 

X de - 4- s (- £ ) [lo (f) + Ii (f ) ] 


sin 3 0 erf (s sin 0) d0 = 


( a+ 50 Il (r)! • 


and the drag coefficient due to specular reflection, referred to the 
frontal projected area, is 

Cd s - — £ «s (- [(3+2s 2 ) Io(|i) + ( l+ 2 s 2 ) li (|i)J (A13) 


Sphere 


As shown in figure 2, the element of area for the sphere is given 


dA = 2itE 2 cos 0 d0 
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Equation (Al) for the incident momentum drag force is 

. P^ 2 , 2.2 cos 8 f .2 2 \ 

dG-i = (nR ) 2s exp(-s sin 81 + 

/7 b 2 l 


Jx sin 8 (2s 2 + 1) erf (s sin 8) 


da 


(A14) 


The integrals are 
if/2 


r /it 

/ ** s cos 8 exp(-s 2 sin 2 S) d8 = — erf(s) 

° / _ 1 \ 

1^/2 y g \ 1 2S 2 ✓ 

/ cos 8. sin 0 e'rf(s sin 8) d8 = ■ 9 x ^ > . S~ s — 1+ — — erf(s) 

J 0 2/ S 2 

and the incident drag coefficient, referred to the frontal projected 
area, is 


CDg - 


202+1 ezp(— a 2 ) + i fr 8 * + 4 ^ .. g -y . erf (s ) 

2s 


vG 3 


‘ (A15) 


s 


The drag force due to molecular re— emission, from equation (A 6), is 

jcR 2 P 2 . Vi! sin 2 8 cos 8 dQ (Al6) 

2 , J 0 s r 

Integration of equation (Al6)' yields the re-emission drag coefficient 

2/TT . 


C;D r 3 s r 

The total drag is given, for diffuse reflection, by 

_ 2 exp(-s 2 ) 1 . 1 1 . . 2</"jt 

C D = — i (1 + — ) + 2(1 + -g. - — r)erf(s) + — — 

Jx s 2s 2 - sr 4s 3s r 


(AIT) 


The drag force. due to speoular reflection on the sphere is given by 
_n2 P ' 4 cos 8 sin 8 


=» Jf/2 

T KR 2 P 

<Jq 






2s sin 8 exp(-s 2 sin 2 8) + 


2 /k ^s 2 sin 2 0 + - ^ erf (s sin 8) 


J 


de 


(ai8) 
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The integrals are 

0 sin £ 0 exp(— e^in 2 ©) d 0 = erf(s) — — - 

- 4s 3 2a 2 

9 ain 0 erf (a sin 0) d0 = SS B ^ fl + erf ( s ) 

2 Vir a \ 4a 2 / 

0 sin 3 0 erf (a sin 0) d0 = 2a - + - \ - exp(-s 2 ) + erf(s) 

8,/its 3 16s 4 

and the drag coefficient due to specular reflection, referred to the 
frontal projected area, is 

°% - Jr? + /—r/) (A19) 


Jt/2 


COS 


l 


5t/2 


COS 


rt /2 


1 


cos 


Cone 


Special attention must-be given to the conical body. Two separate 
cases must be investigated', depending on the relationship of the angle 
of attack of the cone and the semivertex angle of the cone . At an angle 
of attack less than or equal to the semi vertex angle, the entire surface 
area is exposed to the mass velocity, and the total momentum force is 
evaluated by a single integration. For those instances where the angle 
of attack exceeds the semivertex angle, a portion of the surface area is 
shielded from. the mass velocity, and two separate integrations, one for 
the frontal exposed area' and one for the shielded area, are necessary to 
compute the momentum force. 1 Aerodynamic coefficients are then found 
referred to the base area of* the cone. The element of area chosen for 
constant angle of attack as shown in figure 3 is given by 


dA = 


tan 5 h 1 cLcp 
cos S 2 


Case I — Angle of attack less than or equal to the semi vertex 
angle. — The aerodynamic coefficient due to impinging moleoules referred 
to the base area of the cone is given by (from the first term of 
equation (2)) 
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C< = 


* 3/2 s 2 sin 5 


l 


adcp 


(A20) 


where 


a = \^ x PU x +lyPD y +2 z pU z ^) j^exp^p 2 ^ 2 ) +J* pu z [1 + erf (gU x )] 


[1 + erf (pU x )] 


(A21) 


Substitution of the appropriate direction cosines into equation 
(A20) yields the corresponding lift or drag coefficient. 

For drag the direction oosines are 

Z X( i = cos a sin 8 — sin a cos 5 cos cp 


‘'yd 


= — cos a sin 8 — sin a sin 8 cos cp 


> (A22) 


Z zd = sin a sin cp 

and for lift the direction cosines are 

Z X 2 = — sin a sin 8 — cob a cos 8 cos cp 
Zy-j = sin a cos 8 — cos a sin 8 cos cp 
Z Z 2 = cos a sin cp 


} (A23) 


The expressions for the several velocity components are the 
following: 

U x =• U(cos a sin 8 — sin a cos 8 cos cp) 

Uy = -U(cos a 00 s 8 + sin a sin 8 cos cp) 

U z = U sin a sin cp 


(A24) 


At zero angle of attack, equation (A20) reduces to a readily 
integrahle expression and the aerodynamic coefficients for this case 
are 
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s exp(-s 2 sin 2 5) + Vir sin 5 


= 


(il+i) 


1 + erf(s sin 6) 


tfft s 2 sin 5 


(A25) 


and 


Ct„ = 0 


Tlie lift and drag coefficients due to diffuse molecular reflection 
are found, from equation ( 3 ), to he 


71 

^ = 2*ss r sin 5 J l xd xdc P 


(A26) 


and 


2nssv.sin 8,/ X<ic P 


/ 


(A27) 


At zero angle of attack, 


°Dr 


X 


2S S T 


(A28) 


and 

Cl* - 0 

The following relations hold for the total lift and drag coeffi- 
cients for the case of specular reflection: 


Cd S - -T7“ f *xd (PC expC-f 2 ^ 2 ) + 

rf3/ 2 s 2 sin 5 J n l 


4H 


1 + erf(pUx) 


(p^a 


2 1 

x + g 


d<P 


(A29) 


and 
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^ - ,«>. k^i l li1 {** eip( -* ao + 

tfrt 1 + erf(PDx)J + 2 ^ j" ^ 


(A30) 


For zero angle of attack 


„ 2 sin 5 

C-n = — e 


ft 


xp(— s 2 sin 2 5) + 2 (sin 2 5 + — — j [l + erf(s sin 
\ 2s a / 


5)] 


(A31) 


and 


Cls = ° 

Case H — Angle of attack greater than the semi vertex angle.— In 
this case, part of the surface area is shielded from the main flow. 

This shielded area is hounded by straight-line elements of the cone 
and, at these boundaries, the flow is tangent to the surface of the 
cone; that is, the component of Telocity in a direction normal to the 
surface is zero. 

From equation (A2k) 

U x = TJ(cos a sin 5 — sin a cos 5 cos q>x) = 0 

at cp = _ (A32) 


and 


= cos 


tan 5 
tan a 


where <?i denotes the boundary line between the esqposed and shielded 
surface areas and is measured around the cone from the vertical as 
shown in figure k. 


The. total incident momentum for the entire surface area is 



2k 


HACA TN 2423 


Gi = 


pU 2 L a tan S 
2 -fn a a cos & 


< f 1 (l x , PU x »+ly*pTT y »+Z z *pU z *) | esjK-p 2 ^**) 

*/* PU X » [1 - erf(j3U x »)]} lx* [1 - erf ( £U X * ) ] dtp + 

J J(2 x pU x +Z y pU y +Z z pU z ) |exp(-^ 2 U x a ) + Jx pU x [1 + arf(0U x )] j- + 


~ lx [1 + erf(pu x )] 


where 


dcp 


(A33) 


V 

II 

1 

K 

V - - u x 

V 

" hr 

& 

II 

> 

V 

- -*z 

II 

1 

& 

-w.. 


(A34) 


A change In the coordinate axis of the element of area on the 
shielded surface of the cone is the reason for the use of the primed 
quantities. The positive direction of the x axis for the exposed area 
is into the surface; whereas the positive direction of the x axis for 
the shielded area is out of the surface. Since the direction cosines 
have "been derived for the exposed area, reversal of the x axis over the 
shielded area results in the negative quantities. For the previous 
bodies considered, this procedure was unnecessary since both front and 
rear areas were Jointly taken into account. 

The solution of equation ( A33) .» for the incident aerodynamic 
coefficient referred to the base area, yields 



where and ft* are the integrands of equation (A33)» Either the 
lift or drag coefficients can be calculated by substituting the appro- 
priate direction cosines. 


The re-emission drag and lift coefficients for diffuse reflection 


are 


2rtss r s 


(- C 1 **' x ' w + ^” IsaX39 ) 


(A36) 
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and 

CL r - i^k-6 (~j\' ^ 1 ^ XW ) 

The drag coefficient for specular reflection is expressed as 


(A37) 


CDg 


rt 3 / 2 s?sin 5 'f o Zxd ’ { PUxl 

^ [1 - erfOV )] ( W x * 2 


expC-p 2 !^* 2 ) - 
+ |)}d<P + 


^xd -jp'D'x expC-P 2 ^ 2 ) +V it [1 + erfCPU^)] (p 2 ^ 2 + ^j-dcpj (A38) 


The total lift coefficient, Clq, is found by substitution of l x i 
for 2 x d- The values of the above integrals were found by means of 
numerical methods . 
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Figure 2.- Elemental areas for sphere and cylinder : 
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Figure 3.- Element of surface area for cone. 
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Molecular speed ratio, s 



Figure 5- Flat -plate drag coefficient, diffuse reflection. 


Ratio of lift and drag coefficients , 


32 


mCA TN 2423 



Figure 6.- Flat- plate lift-drag ratio, diffuse reflection . 
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Figure 9- Sphere drag coefficient , diffuse and 

specular reflection. 
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Figure 14. - Cone drag coefficient, comparison of free 
molecule and continuum - flow theory ; a* 0°. 
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